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ABSTRACT 

The  subject  of  this  report  is  the  mechanical  behavior  of  stiffened  plates,  basic  structural 
components  of  ships  and  submarines.  The  buckling  loads  of  grillages  subjected  to  axial 
compression  with  and  without  lateral  pressure  are  calculated  using  a  finite  element  based 
eigenvalue  analysis.  Insights  are  obtained  into  the  ways  in  which  the  buckling  loads  and 
modes  vary  with  various  grillage  dimensions. 
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INTRODUCTION 


Stiffened  plates  are  basic  structural  components  of  ships  and  submarines.  These  structures 
are  designed  with  generous  safety  margins  against  overall  collapse  triggered  by  buckling. 
The  object  of  mathematical  modeling  is  to  determine  design  criteria  to  inhibit  buckling  at 
any  stress  less  than  yield.  In  our  earlier  work  (see  References),  we  have  reviewed  the  existing 
literature  and  developed  new  methods  for  predicting  the  buckling  loads  of  simple  plate-beam 
structures.  In  the  present  work,  we  calculate  the  axial  buckling  loads  of  grillages  with  the 
use  of  a  well-known  finite  element  code.  The  goal  is  to  determine  the  effect  on  the  buckling 
loads  of  varying  various  grillage  parameters. 

The  grillage  now  modeled  consists  of  3  base  plates  with  4  longitudinal  and  2  transverse 
stiffeners.  The  dimensions  of  this  grillage  are  denoted  in  the  report  by 

al  =  length  between  transverse  stiffeners 

a2  =  length  between  transverse  stiffener  and  grillage  end 

bl  =  width  between  longitudinal  stiffeners 

b2  =  width  between  outer  longitudinal  stiffener  and  grillage  side  =  bl/3 
tl  =  thickness  of  inner  base  plate 

4 

t2  =  thickness  of  outer  base  plates  =  tl  x  — 

o 

dwl  =  depth  of  longitudinal  webs 
twl  =  thickness  of  longitudinal  webs 
dfl  =  depth  of  longitudinal  flanges 
tfl  =  thickness  of  longitudinal  flanges 
dw2  =  depth  of  transverse  webs 
tw2  =  thickness  of  transverse  webs 
df2  =  depth  of  transverse  flanges 
tf2  =  thickness  of  transverse  flange. 
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The  material  chosen  for  this  study  is  isotropic  steel  with  Young’s  modulus  E  =  3  x  107  psi 
and  Poisson’s  ratio  v  =  .3.  The  imposed  boundary  conditions  are: 

(1)  One  end  of  the  grillage  has  all  3  displacement  components  zero  and  all  3  rotation 
components  zero. 

(2)  The  other  end  of  the  grillage  where  the  force  is  applied  has  axial  displacement 
constant  with  the  other  2  displacement  components  zero  and  all  3  rotation  com¬ 
ponents  zero. 

(3)  The  ends  of  the  transverse  stiffeners  have  vertical  displacement  zero  but  the  other  2 
displacement  components  and  all  3  rotation  components  are  free. 

(4)  All  other  nodes  on  the  plates  edges  are  completely  free. 

Because  of  the  boundary  conditions  and  the  fact  that  the  thickness  t2  of  the  outer  base  plates 
is  chosen  to  be  |  times  the  thickness  tl  of  the  inner  base  plate,  the  buckling  deflections  of 
the  outer  bays  are  usually  small  compared  to  the  buckling  deflections  of  the  central  bay. 
Thus  we  calculate  the  buckling  stress  in  the  central  bay,  defined  as  the  initial  end  force 
divided  by  the  cross-sectional  area  of  the  central  bay,  even  though  actual  buckled  members 
may  be  elsewhere. 

The  finite  element  code  used  is  MSC  Nastran/Patran.  The  pre-processing  and  post¬ 
processing  is  done  with  Patran  (versions  1.4  and  1.5).  The  buckling  analysis  (solution  105) 
is  done  with  Nastran  (version  68).  The  base  plates,  flanges,  and  webs  are  each  modeled  with 
Quad  4  plate  elements.  The  mesh  size  is  denoted  in  this  report  by 

ml  =  mesh  length  of  all  elements,  mesh  width  of  base  plate  elements 
m2  =  mesh  width  of  web  and  flange  elements. 

That  is,  the  base  plating  has  square  elements  of  size  ml  x  ml,  while  the  webs  and  flanges 
have  rectangular  elements  of  size  ml  x  m2.  Note  that  this  is  a  full  finite  element  model  of 
the  entire  grillage  and  does  not  assume  symmetry.  We  use  linear  geometry  for  the  axial  load 
only  cases,  nonlinear  geometry  for  the  pressure  cases,  and  linear  material  properties  for  all 
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cases.  All  computations  are  performed  on  a  Silicon  Graphics  Indy  Workstation. 

The  first  phase  of  our  work  is  to  compare  the  finite  element  model  with  experimental 
results  on  a  small  scale  grillage  under  axial  and  pressure  loading.  The  second  phase  is  to 
vary  the  dimensions  on  a  full  scale  model  to  assess  changes  in  buckling  stress  and  mode 
under  axial  load  only. 
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EXPERIMENTALLY  TESTED  GRILLAGES 

We  first  consider  the  grillages  (  Figure  1)  tested  in  the  USNA  Ship  Structures  Laboratory 
testing  rig  (Figure  2).  We  choose  the  following  values  for  the  model  parameters: 


al  =  36" 
a2  =  36" 
bl  =  9" 


dwl  =  3" 


twl  =  .114" 
dfl  =  1.844" 
tfl  =  .171" 
dw2  =  9" 
tw2  =  .25" 
df2  =  3" 
t£2  =  .25" 
ml  =  1" 
m2  =  .5". 

We  find  that  the  buckling  load  is  not  significantly  changed  by  adopting  smaller  values  for 
the  mesh  sizes  ml  or  m2. 

The  prebuckling  stage  is  approximately  one  of  uniform  axial  compression  with  little 
bending  (Figure  3).  The  grillage  buckles  under  a  force  of  470  kips.  The  buckling  mode  is 
a  deformation  involving  primarily  bending  with  little  stretching  (Figure  4).  The  inner  base 
plate  buckles  into  a  square  quilt  of  half- wavelength  9"  (Figure  5).  The  webs  (both  transverse 
and  longitudinal)  also  buckle  with  half- wavelength  9"  (Figures  6-7).  The  flanges  mainly  just 
twist  about  their  center  line.  We  call  this  buckling  mode  TRIPPING  to  be  consistent  with 
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our  earlier  papers.  * 

We  also  ran  nonlinear  analyses  (solution  106)  for  the  cases  of  bottom  pressure  p  =  0, 
p  =  5  psi,  or  p  =  20  psi.  For  the  case  p  =  0,  the  plate  only  compresses  in-plane  and  does  not 
buckle.  However,  the  buckles  are  clearly  visible  by  the  time  the  axial  load  reaches  475  kips 
for  the  p  =  5  psi  case  and  450  kips  for  the  p  =  20  psi  case.  But  since  the  buckling  pattern 
develops  gradually  as  the  axial  load  increases,  it  is  hard  to  pin  down  an  exact  “buckling” 
load.  It  is  interesting  that  the  middle  plate  develops  only  3  half  waves  (instead  of  4  as  in 
the  buckling  analysis  for  p  =  0).  The  magnitude,  of  displacement  for  the  p  =  20  psi  case  is 
shown  in  Figure  8  with  no  axial  load,  and  in  Figure  9  with  550  kips  axial  load.  If  we  increase 
the  load  much  further  than  500  -  550  kips,  the  displacements  become  large  and  the  code  will 
not  give  a  solution  (Figure  10). 

A  table  comparing  model  results  to  experimentally  measured  collapse  loads  for  6  nomi¬ 
nally  identical  grillages  is  given  below: 


Model  (kips) 

Experiments  (kips) 

0 

470 

326,  312,  and  301 

5 

less  than  475 

316 

10 

not  calculated 

306 

20 

less  than  450 

296 

*In  our  earlier  work,  we  did  not  allow  the  stiffener  cross-section  to  deform  in  its  plane.  This  assumption 
led  to  a  predicted  buckling  load  which  was  much  too  high. 
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Note  that  both  theory  and  experiments  indicate  that  the  presence  of  lateral  pressure 
(on  the  plating  side)  tends  to  weaken  the  grillages.  The  disparity  between  the  buckling 
loads  of  the  perfect  model  and  the  experimental  measurements  is  probably  due  to  geometric 
imperfections,  residual  stress,  differing  boundary  conditions,  nonlinear  material  properties, 
etc.  in  the  tested  grillages. 
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Hydraulic 

Cylinder 


USNA  Ship  Structutes  Laboratory  testing  rig. 


FRINGE  PLOT  LC=1 .23  RES=2.1(VEC-MAG)  MSC/PATRAN  R-1 .4  MSC/NASTRAN  29-Jan-96  1 1 :55:35 
DEFORMATION  PLOT  LC=1 .23  RES=2.1  MSC/PATRAN  R-1 .4  MSC/NASTRAN  29-Jan-96  1 1 :55:35 

Prebuckling  magnitude  of  displacements 
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FRINGE  PLOT  LC=2.24  RES=3.1(Z-COMP)  MSC/PATRAN  R-1.4  MSC/NASTRAN  29-Jan-96  13:31:39 
DEFORMATION  PLOT  LC=2.24  RES=3.1  MSC/PATRAN  R-1 .4  MSC/NASTRAN  29-Jan-96  13:31 :39 
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COMP)  MSC/PATRAN  R-1 .4  MSC/NASTRAN  29-Jan-96  13:40:05 
>=3.1  MSC/PATRAN  R-1 .4  MSC/NASTRAN  29-Jan-96  13:40:05 
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FRINGE  PLOT  LC=3.90  RES=2.1(VEC-MAG)  MSC/PATRAN  R-1.4  MSC/NASTRAN  12- Feb-96  10:37:53 
DEFORMATION  PLOT  LC=3.90  RES=2.1  MSC/PATRAN  R-1 .4  MSC/NASTRAN  12-Feb-96  10:37:53 
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Figure  10 


GRILLAGES  OF  VARIOUS  DIMENSIONS 


Generic  model 

We  define  the  generic  model  as  a  grillage  with  the  following  parameter  values: 


al  =  96" 

a2  =  48" 

bl  =  24" 


dwl  =  6" 

twl  =  .2" 

dfl  =  4" 
tfl  =  .225" 

dw2  =  12" 

tw2  =  .2" 

df2  =  4" 

t£2  =  .225" 

ml  =  2" 

m2  =  1". 
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The  buckling  stress  in  the  central  bay  of  the  generic  model  is  66.2  ksi.  The  buckling 
mode  is  TRIPPING  of  the  central  bay  (Figure  11).  The  half- wavelength  of  the  buckles  is 
24". 

Note  that  the  outer  bays  of  the  generic  model  are  made  half  as  long  as  the  inner  bay,  to 
reduce  the  size  of  the  finite  element  model.  The  buckling  load  of  the  full  scale  model  differs 
from  that  of  the  generic  model  by  an  insignificant  amount. 

For  the  remainder  of  the  report,  we  consider  various  grillages  with  only  one  dimension 
different  from  those  of  the  generic  model. 

Variation  of  al,  length  between  transverse  stiffeners 

We  consider  6  grillages  having  differing  values  of  al.  All  of  the  other  parameters  are  the 
same  as  for  the  generic  grillage,  except  for  the  case  al  =  144"  we  choose  the  mesh  length 
ml  =  4"  to  make  the  program  size  smaller.  The  buckling  stress  decreases  monotonically 
with  al  (Figure  12).  The  buckling  mode  is  TRIPPING  in  all  cases.  For  the  5  cases  with 
al  a  multiple  of  24",  the  half-wavelength  of  the  buckles  is  24".  For  the  case  al  =  60",  the 
half- wavelength  is  30"  (Figure  13). 

Variation  of  bl,  width  between  longitudinal  stiffeners 

We  consider  4  grillages  having  bl  differing  multiples  of  12".  The  buckling  stress  decreases 
monotonically  with  bl  (Figure  14).  The  buckling  stress  is  somewhat  larger  than  the  classical 
buckling  stress  of  a  square  simply-supported  plate  (Figure  15): 

£7T2(tl)2 
3(1  —  j/2)(61)2 

The  buckling  mode  is  TRIPPING  (similar  to  Figure  11)  and  the  half-wavelengths  of  the 
buckles  is  bl. 
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Variation  of  tl,  thickness  of  inner  base  plate 

We  consider  5  grillages  having  differing  values  of  tl.  The  buckling  stress  reaches  a  maximum 
value  for  an  optimal  tl  (Figures  16  -  17).  For  the  3  lowest  cases,  the  buckling  mode  is 
TRIPPING,  with  the  half- wavelengths  being  19.2"  for  the  case  tl  =  and  24"  for  the 
cases  tl  =  |"  and  tl  =  jg"  For  the  cases  tl  =  |"  and  tl  =  j|"  the  buckling  mode  is 
GLOBAL.  The  entire  central  bay  buckles  in  (Figure  18),  while  the  stiffeners  deflect  inwards 
and  laterally  (Figure  19). 

Variation  of  dwl,  depth  of  longitudinal  webs 

We  consider  7  grillages  having  differing  values  of  dwl.  The  buckling  stress  reaches  a  max¬ 
imum  value  for  an  optimal  dwl  (Figure  20).  For  the  cases  dwl  =  3"  and  dwl  =  4",  the 
buckling  mode  is  GLOBAL  (Figure  21).  For  the  cases  dwl  =  5",  dwl  =  6",  and  dwl  = 
9",  the  buckling  mode  is  TRIPPING  with  half-wavelength  of  24".  For  the  cases  dwl  =  10" 
and  dwl  =  11",  the  buckling  mode  is  LOCAL  (Figure  22).  The  web  and  flange  of  a  middle 
stiffener  of  the  central  bay  buckle  with  a  half- wavelength  dwl.  We  find  that  the  critical  load 
for  the  local  buckling  case  is  not  significantly  changed  by  adopting  smaller  values  for  the 
web  and  flange  mesh  element  width  m2. 

Variation  of  twl,  thickness  of  longitudinal  webs 

We  consider  7  grillages  having  differing  values  of  twl.  The  buckling  stress  increases  mono- 
tonically  with  twl  (Figure  23).  For  the  case  twl  =  .05",  the  web  of  a  middle  stiffener  of  the 
forward  bay  buckles  LOCALLY,  and  for  the  case  twl  =  .1",  the  web  of  a  middle  stiffener  of 
the  central  bay  buckles  LOCALLY  (Figure  24).  For  the  cases  twl  =  .15",  .2",  .3",  and  .7" 
the  buckling  mode  is  TRIPPING  with  half-wavelength  24"  (Figure  25).  For  the  case  twl  = 
1",  the  buckling  mode  is  TRIPPING  with  half-wavelength  19.2". 


Variation  of  dfl,  depth  of  longitudinal  flanges 

We  consider  4  grillages  having  differing  values  of  dfl.  For  the  cases  dfl  =  2",  4",  and  6", 
the  buckling  mode  is  TRIPPING  and  the  buckling  stress  varies  little  (Figure  26).  For  the 
case  dfl  =  0  (no  flange),  the  web  of  a  middle  stiffener  of  the  central  bay  buckles  LOCALLY 
(Figure  27). 

Variation  of  tfl,  thickness  of  longitudinal  flanges 

We  consider  4  grillages  having  different  values  of  tfl.  For  the  cases  tfl  =  .225",  .5",  and  1", 
the  buckling  mode  is  TRIPPING  and  the  buckling  stress  varies  little  (Figure  28).  For.  the 
case  tfl  =.l",  the  flange  of  a  middle  stiffener  of  the  central  bay  buckles  LOCALLY  (Figure 
29). 

Variation  of  dw2,  depth  of  transverse  webs 

We  consider  3  grillages  having  differing  values  of  dw2.  For  all  3  cases  dw2  =  4",  12",  and 
16",  the  buckling  mode  is  TRIPPING  and  the  buckling  stress  varies  little  (Figure  30). 

Variation  of  tw2,  thickness  of  transverse  webs 

We  consider  3  grillages  having  differing  values  of  tw2.  For  the  cases  tw2  =  .15"  and  .2",  the 
buckling  mode  is  TRIPPING  and  the  buckling  stress  varies  little  (Figure  31).  For  the  case 
tw2  =  .1",  a  web  buckles  LOCALLY  under  a  tensile  force  (Figure  32).  Compression  in  the 
web  arises  from  Poisson  ratio  effects. 

Variation  of  df2,  depth  of  transverse  flanges 

We  consider  2  grillages  having  differing  values  of  df2.  For  both  cases  df2  =  0  (no  flange) 
and  df2  =4",  the  buckling  mode  is  TRIPPING  and  the  buckling  stress  varies  little  (Figure 
33).  Note  that  eliminating  the  transverse  flange  would  not  reduce  the  buckling  strength  of 
the  grillage  under  axial  compression. 
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Variation  of  t£2,  thickness  of  transverse  flanges 

We  consider  2  grillages  having  differing  values  of  tf2.  For  both  cases  tf2  =  .1"  and  tf2  =.225", 
the  buckling  mode  is  TRIPPING  and  the  buckling  stress  varies  little  (Figure  34). 

Summary 

Through  study  of  all  these  different  results,  we  can  gain  insights  into  the  ways  in  which  the 
buckling  loads  and  modes  vary  with  the  various  dimensions. 

The  buckling  stress  decreases  less  with  length  al  than  with  width  bl  between  stiffeners. 
The  buckling  stress  is  larger  for  an  optimal  value  of  plating  thickness  tl  and  depth  dwl  of 
longitudinal  webs.  The  buckling  stress  increases  with  thickness  twl  of  longitudinal  webs. 
The  buckling  stress  is  not  much  affected  by  increases  in  the  dimension  dfl  and  tfl  of  the 
longitudinal  flanges,  dimensions  dw2  and  tw2  of  transverse  webs,  or  dimensions  df2  and  tf2 
of  transverse  flanges.  However,  the  buckling  stress  is  smaller  for  too  small  values  of  thickness 
twl  of  longitudinal  webs,  depth  dfl  and  thickness  tfl  of  longitudinal  flanges,  and  thickness 
tw2  of  transverse  webs. 

We  have  found  three  distinct  failure  modes  which  we  have  named  TRIPPING,  GLOBAL, 
and  LOCAL.  For  the  majority  of  our  cases,  the  failure  mode  is  TRIPPING  with  the  inner 
base  plate  buckling  into  a  square  quilt,  which  forces  the  webs  to  buckle  into  a  rectangular 
quilt  with  the  same  wavelength,  while  the  flanges  twist  (Figure  11).  However,  if  a  square 
pattern  cannot  fit  into  the  base  plate  (Figure  13),  or  if  the  relative  strength  of  the  stiffeners 
is  much  greater  than  the  base  plate  (Figure  25),  the  inner  base  plate  may  buckle  into  a  non¬ 
square  quilt,  and  the  flanges  may  also  undergo  lateral  deformation.  If  the  relative  strength 
of  the  inner  base  plate  is  much  greater  than  the  stiffeners,  the  failure  mode  may  be  GLOBAL 
with  the  entire  central  bay  buckling  in  (Figures  18-19  and  21).  If  the  relative  strength  of  a 
stiffener  is  too  small,  it  may  collapse  LOCALLY  before  any  of  the  other  components  buckle 
(Figures  22,  24,  27,  29,  and  32). 
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Figure  34 


CONCLUSIONS 


We  can  think  of  the  buckling  stress  as  being  represented  by  a  surface  in  n-dimensional 
parameter  space  (n  =  11  in  our  studies).  The  surface  is  continuous  but  has  discontinuous 
slope  at  the  boundaries  between  subsurfaces  corresponding  to  TRIPPING,  GLOBAL,  or 
LOCAL  modes.  Our  curves  of  buckling  stress  versus  a  parameter  are  the  intersections  of  the 
surface  with  vertical  planes  passing  through  a  given  point  on  the  surface. 

To  fit  the  entire  surface  with  an  analytical  approximation  would  appear  to  be  a  difficult 
computational  problem  requiring  the  calculation  of  numerous  points  on  the  surface.  The 
parameter  dependence  of  the  buckling  load  is  coupled,  i.e.,  all  the  parameters  would  have  to 
be  varied  simultaneously. 

In  future  work  the  finite  element  eigenvalue  code  could  be  used  to  design  the  ideal  struc¬ 
ture.  One  possible  design  criterion  would  be  to  find  the  grillage  of  minimum  weight  among 
all  possible  grillages  with  buckling  stress  equal  to  the  yield  stress.  The  candidate  points 
would  lie  on  the  intersection  of  the  surface  with  a  horizontal  plane. 
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